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We show that the surface states of pristine 3D topological insulators (TI) are analogs of ferromag¬ 
netic half metals due to complete polarization of an emergent momentum independent pseudospin 
(SU(2)) degree of freedom on the surface. To put this claim on firm footing, we present results for 
TI surfaces perpendicular to the crystal growth axis, which clearly show that the tunneling con¬ 
ductance between two such TI surfaces of the same TI material is dominated by this half metallic 
behavior leading to physics reminiscent of a spin-valve. Further using the generalized tunnel magne- 
toresistance derived in this work we also study the tunneling current between arbitrary TI surfaces. 

We also perform a comprehensive study of the effect of all possible surface potentials allowed by 
time reversal symmetry on this spin-valve effect and show that it is robust against most of such 
potentials. 


I. INTRODUCTION 

Three dimensional topological insulators (3D TI)^“^ represent a distinct class of 3D band insulators which have 
topologically protected metallic surface states. Nontrivial topology of band structures of these materials leads to 
the immunity of these surface states against a variety of disorder and interaction potentials. Possibility of surface 
magnetization and related spin textures of Fermi surface in these materials due to external doping^“^^ has also been a 
topic of great interest. Here we show that, even the pristine TI surface can act like a ferromagnetic half metad^, though 
not due to spin but due to an emergent pseudospin degree of freedom^^. A natural way to confirm a half-metallic 
behavior is an observation of perfect spin orthogonality (complete suppression of tunneling due to orthogonality of 
the states in the spin sector) in tunnel conductance between two half-metals with opposite directions of polarization 
and then a lifting of the orthogonality by tilting the polarization direction of one half-metal with respect to another. 
The central result of this work is the discovery of such (pseudo) spin orthogonality physics and its subsequent lifting 
due to the rotation of pseudospin polarizations on the surface states of 3D TIs. 

The low energy physics of popular TI materials like Bi 2 Se 3 can be described by a four band model arising from 
two SU(2) degrees of freedom.Of these two SU(2) degree of freedom, one of them is dispersing^^ while the 
other one is non-dispersing. This fact is natural as the surface state itself comes into being by freezing (making it 
momentum independent) two out of four bulk degrees of freedom (four band model) on the planar boundary. This 
makes the frozen SU(2) degree of freedom fully polarized and hence surface state is endowed with a half-metal like 
character. 

It is interesting to note that the polarization of the frozen degree of freedom is exactly opposite on the opposite 
surfaces for any angle that the surfaces make with the crystal growth axis, hence leading to distinction between 

the top and bottom surfaces which is very similar to the distinction of north and south poles of a bar magnet. It can 
be understood from the observation that, when a bar magnet is broken into two halves, it exposes two new ends with 
opposite polarity such that each new bar magnet has opposite polarity at its two ends. In an identical fashion when a 
sample of TI with exposed surfaces corresponding to 0 and OPir is sliced parallel to the plane creating an angle 0 with 
the crystal growth axis, creating two new TI films, the newly exposed surfaces will have exactly opposite polarizations 
in the frozen degree of freedom consistent with the scenario of a bar magnet. The fact that the different surfaces 
of the TI are very distinct from one another due to this frozen degree of freedom has been mostly underestimated 
and ignored. One of the most important outcome of our analysis is, though frozen, these degree of freedom can have 
strong influence on transport across junctions of TI via physics analogous to the spin-valve effect. In our analysis, the 
role of the edge states which may appear at the edges of a 3D sample^^’^^ will have a negligible effect and hence are 
neglected. 

The rest of the paper is organized is follows, in Sec. II we describe the model Hamiltonian for the surface states 
their pseudospin textures. In Sec. HI we sketch the calculation of the tunneling current between two TI surfaces 
and derive the expression for the generalized tunnel magnetoresistance followed by a discussion of identification of 
the half-metallic behavior via tunneling current in Sec. IV. In Sec.V we discuss the effect of interplay of the general 
tunnel magnetoresistance and the shape of Fermi surface on tunneling current between two arbitrary surfaces. The 
robustness of the half-metallic behavior against surface potentials is discussed in Sec. VI followed by concluding remarks 
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FIG. 1. (a) and (b): A schematic showing the two different junctions A and B of TI surfaces with an insulating barrier between 
them. The spin textures [(c) and (d)] and the orbital pseudospin textures [(e) and (f)] of the two surfaces are plotted on the 
Dirac cones of the surface states. The orbital pseudospin textures are completely polarized along the positive (blue) or negative 
(red) X direction of the r-space. 


in Sec.VIL 


II. MODEL AND PRELIMINARIES 

In this section, we describe the model Hamiltonian for the surface states of the 3D TI. We sketch the derivations of 
the surface state spinors and their properties, setting up consistent notations used throughout the rest of the paper. 
For details refer to Refs. [5, 13, 14, and 17]. 

Popular 3D TI material crystals like Bi 2 Se 3 have a quintuple layered unit cell comprising of two equivalent Se layers 
(Sel and Sel'), two equivalent Bi layers (Bil and Bil') and one inequivalent Se layer (Se2). Both, Bi and Se have the 
p-orbitals as their outermost shells, hence it is natural to consider only the p-orbitals for the electronic properties. As 
the hybridization of the Bi and Se orbitals, formation of the anti-bonding and bonding orbitals, crystal field splitting, 
and the effect of spin-orbit coupling are taken into account, it turns out that there are four states close to the Fermi 
level which participate in the low-energy description of the materials. We briefly discuss the effect of each of these 
steps in the following: 

• The hybridization of the Bi and Se orbitals pushes up the Bi energy level forming new hybridized degenrate 
levels denoted by B and B' , where as the Se energy levels are pushed down forming new hybridized levels 
denoted as S', S' and SO, out of which S and S' are degenerate. 

• Among the four states near the Fermi level two of them arise from the bonding orbital between the B and H', 
denoted as P1+, and the other two states arise from the anti-bonding orbital between S and S' denoted as P2~. 
Note that due to inversion symmetry, the bonding and anti-bonding orbitals have definite parity. 

• Each of the two states coming from P1+ or P2~ have equal and opposite total angular momentum (sum of 

orbital angular momentum and spin) which turns out to be ±1/2. The states with total angular momentum 
±1/2 are constituted by linear combination of Pz orbitals with t spin and ± ipy orbitals with | spin. 

Similarly a linear combination of Pz orbitals with | spin and p- = Px — ipy orbitals with t spin form the states 
with total angular momentum —1/2. 

• Hence the four states near the Fermi energy can be represented as |Pl+,±^) and |P2“,±^). Note that this 
naturally provides us with two SU(2) degrees of freedom to represent the states. The first SU(2) degree of 
freedom represented by P1+ and P2~ have different orbital parities, hence we call it the orbital pseudospin and 
denote it by the set of Pauli matrices r. It also turns out that for the surface states, the spin angular momentum 
is proportional to the total angular momentum^^, hence we call the other SU(2) degree of freedom, the total 
angular momentum simply the spin and denote it by the set of Pauli matrices cr. 
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As mentioned in Sec. I, the surface states are labeled by an angle 0 between the crystal growth axis and the normal to 
the surface. For the simple case of 6> = 0 (top surface) and 0 = it (bottom surface), the dispersing and frozen degree of 
freedom are given by cr and r operators respectively. The direction of spin, cr, is determined via the spin-momentum 
locking angle while the r points solely along the positive x (negative x) direction for the top (bottom) surface^^’^^ 
for all momentum states hence defining a r polarized half metal. 


A. Surface states 


In a coordinate system, where the exposed TI surfaces lie parallel to the x-y plane with its crystal growth axis at 
an angle 0 to the ^-axis, the composite Hamiltonian governing both the opposite surfaces can be written in the form 

U,{d)=viikyS^-vik:,Sl ( 1 ) 

where S is the dispersing SU(2) degree of freedom given by^^’^^ 


So = - ^OT^ ^ + ^oT^ ^ (2) 

with ao = Vz cosO/vs, Po = s'mO/vs, {v^ cos 0)^ + sin^)^, and cos^Tcr^^^^ sin^. The above 

surface Hamiltonian features a pair of degenerate Dirac-like spectra given by the dispersion 


Ek,±{^) = E^Jvjk^+vf^kl . 


( 3 ) 


Note that the two Dirac cones correspond to the surface states on the two opposite surfaces, each surface has only 
one Dirac one, as these are strong TIs. 

In the T 0 a basis Eq.(l) assumes a block diagonal form. 


ns{0) 


fcT.Btie) 0 \ 

1 , 0 <T.B-{e)) ’ 


where the effective magnetic fields are given by 

= {kyV\\ {a cos sin^), —k^vi^ 

kyV\\ (=F/d cos 0 ^ a sin 0}. 


Parametrizing these effective magnetic fields as 

■Sfc = l-Bfc |{sin6>± cos</>±,sin6>± sin</>±,cos6>±}, 


( 4 ) 

( 5 ) 

(6) 


the two eigenstates coming from each of the two blocks can be written as 


Ixl+W) = 


cos 


; lx|-(fe))=( y 

\cos ^ 


( 7 ) 


where |x^^) and \Xe-) are the positive and negative energy eigenstates respectively of the first block, and and 

Ix^-) are of the second block of Eq.(4). The eigenstates of the Hamiltonian in Eq.(4) can then be trivially expressed 
as 



( 8 ) 


Note that the two-fold degeneracy of the spectra in Eq. 3 comes due to the fact that = |^^|. 

The correct spinors describing the surface states on the two surfaces can be constructed by taking appropriate 
linear combinations of the states of the same energy. Obviously, the coefficients in the linear combinations can differ 
only by phase which we denote by (j)R. This phase is fixed by demanding the state to be an eigenstate of 
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FIG. 2. The orbital texture for three different surfaces corresponding to ^ = 0,7r/4, and7r/2 are plotted in the momentum space. 
The axes of the figures are kx and ky which are the local in-plane axes. It is important to note that for the ^ = 0 surface, 
the orbital pseudospin shows a complete polarization in the pseudo x direction, however as one moves towards a more oblique 
surface, this polarization begins to go down and at the same time, the polarization in the pseudo y direction begins to pick up 
the texture. 


with eigenvalue ±1 (for the surface state of 0 and tt + ^ respectively) where represents the pseudospin degree of 
freedom given by 


To = {aoT^ + l3eTy G) (8) r^}. 


( 9 ) 


Hence the very construction for surface state for any 0 naturally involves imposing a complete polarization on all 
the surface electron states independent of their momentum, making it a perfect analog of an ferromagnetic half metal 
where spins are fully polarized independent of momentum. We refer to these states respectively as the top and the 
bottom surface corresponding to a particular 0. This leads to the expression of (pR as 


= +(_) 


The resultant surface states are written as 


cos(6>+/2) 


acos{0~/2) — Psm{0~/2)e^^ 


t(b) 


( 10 ) 


( 11 ) 


which ultimately leads to the surface states as 

where t(b) in the superscript denotes the top(bottom) and c{v) in the subscript denotes the conduction(valence) band. 


B. Pseudospin textures 

The spin and orbital texture on the surface can be measured by just sandwiching the respective operators between 
the correct surface state spinors and calculating the expectation value. The spin expectation is given by (cr) = 
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{{Xe+W\Xe+) + {Xe+\^\Xe +))following expressions for the individual components 


{a^){k) = 

(cr^)(fe) = 

{a^){k) = 


kyV\\a cos 9 
^v'^kl+vlkl 
-kxVi 

^vjkl+vlkl 

k.jVwasinO 


(13) 


^vjkl+vlkl 

By a similar calculation, the orbital pseudo spin texture in the conduction band can be also be formally expressed as 

(r) = Re[e*'*^^Ti2(x^+|xE+)] (14) 

Using the form of the eigenstates the component-wise expectation values can be obtained as^^ 

Vz COS 0 


{tx) = cos0i?sin6>fc = 


^3 


{Ty) = sin(/)i^ sin6>fc = 


—VzV?\kr sin^ 




-vm 


(15a) 


(15b) 


{tz) = 0. (15c) 

A few interesting features in the orbital pseudospin texture are of note. Firstly, as in the case of the spin texture, 
the magnitude does not stay constant over the Fermi surface, however, the expectation of (r^,) stays the same over 
the momentum space, though the value it takes depends on the surface taken. On the other hand, for any arbitrary 
surface not perpendicular to the crystal growth axis, (r^) has a texture. Secondly, it is also intersting to note that the 
texture of the magnitude of the pseudospin in the momentum space mirrors that of the spin texture. Thirdly, the {jz) 
is identically 0. This is a manifestation of the fact that on any surface, the contribution to the state coming from the 
two parity orbital sectors can differ only by a phase. This can also be seen from the form of the surface state spinor 
(Eq.(12)). Indeed in a low energy field theoretical description of the surface, this is what should be expected. These 
features are shown in Fig. (2). As in Fig. (2), the first row shows the surface which has the trivial orbital texture with 
the orbital pseudospin pointing unifomly in the x direction in the pseudospin space. 

The relative textures of two TI surfaces play a central role in the generalized magnetoresistance to tunneling current 
between them as we show in the next section. 


III. GENERALIZED TUNNEL MAGNETORESISTANCE 

In this section we sketch the calculation of the tunneling current between two TI surfaces from which the form of the 
generalized tunnel magnetoresistance naturally comes out. We consider a translation invariant tunnel Hamiltonian 
describing tunneling between two planar surfaces of TI with arbitrary 0. The unperturbed Hamiltonian of the 
combined system (of the two TFs) can be written as 

Ho = y dk (-E'i(k, 0i)cj + H2(k, 6*2)4,k®2.k), (16) 

where Ci/ 2 ,k is the fermionic annihilation operator for the two TI surfaces (with 0i and O 2 ) and expression for F' 1/2 is 
given in Eq.(3). Translation invariance leads to a momentum resolved tunnel-Hamiltonian density given by 

Ht = J dk {^k4,kC2,k + 2;k4.k®i,k}> (17) 


where z\^ is the overlap between the two surface state wave functions at momentum k. 
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The resultant current at time t can be written as 

{i{t)) = {GS{t)\i{t)\GS{t)) 

= {9\f dt' (18) 

J — oo 

where \GS{t)) and \g) are respectively the time-dependent and time-independent ground states of the combined system 
and the last line in the above equation is obtained keeping the leading term in J in the Dyson series of the evolution 
operator^^. Using the algebra of the fermionic operators and time translation invariance we arrive at 

/2p poo n 

{i{t oo)) := (I) = —{gl J dt j dk|2;kp[6],(i)6k(0)ak(i)ak(0) “ ^k(0)6],(i)ai.(0)ak(i)] |fl), (19) 

which finally leads to the following expression of the tunneling current density given by 

{I) = ^ Idk |%|"AnF(k, II,, 112 ) X jdE[5iE,ik) - E)5iE2ik) - E)]. (20) 

Here AnF(k,/ii,/i 2 ) = ^F(^i(k),/ii) — nF(^ 2 (k),/i 2 ) is the difference in the Fermi functions of the two surfaces and 

= [1 + (<^)i • (<^)2 + (t)i • (t) 2 + 0(7j)i(ri 0crj)2]/4. (21) 

hi 

Eq.21 is the expression of the generalized tunnel magnetoresistance for tunneling between two TI surfaces and is one 
of the main results of the work. This decomposition is both nontrivial and instructive as it allows for interpretation of 
the total tunneling current in terms of individual responses from cr, r sectors and a correlation term between them. 
The surface states being eigenstates of Tq = agr^ -h (g) lead to a momentum independent part in jzkp given 
by \zf\‘^ = {rx)i{rx )2 ^ {Ty 0 (Ty)i{Ty 0 ay) 2 ^ which is representative of the overlap between the frozen sectors of the 
respective surfaces which does not depend on the momentum. These frozen sectors are crucial towards identifying 
the half-metallic degree of freedom on the TI surfaces as we shown in the next section. 


IV. IDENTIFYING THE HALF-METALLIC DEGREE OF FREEDOM 

In this section we explicitly identify the half-metallic degree of freedom via tunneling current measurements. We 
specifically choose to work with junctions of ^ = 0 and tt surfaces, as and reduce to the spin a and the orbital 
pseudospin r respectively. Consequently the expression for the general tunnel magnetoresistance (21) simplifies to 

^ |^ l + (^)i(k)-H2(k) ^ j^ l + (rh-(r)2 ^ ^ (22) 

The first term in Eq.22 represents the momentum dependent overlap of spin orientations of the spin-momentum 
locked states on the two surfaces. The second term is of more interest here as it encodes the momentum-independent 
overlap of the pseudospin sectors of the surface states and the robust half-metallic beahvior can arise from this sector. 
Now, we know that r points along positive and negative x-axis for = 0 and 0 = tt surfaces for a given TI crystal. 
Hence for a junction between ^ = 0 surface of one TI with the 0 = tt surface of another TI, the first term in Eq.(22) 
vanishes identically ((t)i • (t )2 = — 1 ) leading to a perfect orthogonality in the pseudospin sector of the wavefunctions 
which for the ^ = 0, tt surfaces is the pure orbital pseudospin sector. This directly stems from the half metallic 
behavior discussed above (see Junction-A in Eig.(l)). 

Next we note that, if both the participating surfaces at the junction are ^ = 0 surfaces for their respective crystals 
(see Junction-B in Eig.(l)), the orbital orthogonality is lifted ((t)i • (t )2 = 1). But, even if we start with two identical 
crystals with their crystal growth axis being parallel (taken as z-axis) to each other, still we have to rotate the crystal 
growth axis of one of the crystals with respect to the other by tt so that the ^ = 0 surfaces of the two crystals could 
face each other hence forming the junction. A careful observation immediately reveals that the tt rotation of the 
crystal growth axis leads to a perfect orthogonality between the spins (cr) of the electron states on the two surfaces 
for a given momentum. Hence, though the orbital orthogonality gets lifted in the Junction-B, an orthogonality spin 
sector ((cr)i(k) • (cr) 2 (k) = —1) arises leading to zero tunneling current. So, we have zero tunneling current for 
both Junction-A and -B. This fact makes the possibility of distinguishing the parallel (Junction-B) from anti-parallel 
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FIG. 3. Left: The color density plot shows the out of plane magnetization ((cr^)) and the vectors show the in-plane x and 
y components for the surface states in the presence of warping. Right: The figure shows that how the spin orthogonality 
stays intact even after the orbital pseudo spin orthogonality is lifted by rotating the sample by tt such that both the surfaces 
participating in the junction have their orbital pseudo spin pointing uniformly in the same direction. The red and blue shades 
in the surface refer to the states having (r^) = ±1 respectively and the vectors show the spin texture on the Fermi surface. 


(Junction-A) orbital pseudospin configuration obscure as the absence of the orbital orthogonality in Junction-B is 
masked by the orthogonality in the spin sector. 

The half-metallic behavior manifests itself as the orbital orthogonality in Junction-A, hence to contrast its presence 
and absence in Junction-A and -B respectively, the spin {a) orthogonality has to be lifted. This can be achieved 
by strong doping which is large enough to push the Fermi surface into the hexagonally warped regime of the energy 
spectrum^^. Next, we demonstrate how this can be used to lift the spin orthogonality. 

The low energy Hamiltonian for the TI surfaces to next higher order after linear is a term cubic in momentum 
which couples to the out-of-plane magnetization. The Hamiltonian with the inclusion of the warping term is given by 

■Hw = (cr X k)^ + + fc^)(T^. (23) 

Such warping term induces a 3-fold symmetric pattern in the out-of-plane spin polarization (cr^), with the Fermi 
surface breaking up into six symmetric regions of alternating positive and negative (cr^) as shown in the left panel of 
Fig. (3). The momenta with polar angle 6 >k* = ^ 7 r/ 3 , where n = 0,1... 5, have maximal (a^) with its sign given by 
(-!)”• 

The presence of the out-of-plane magnetization allows us to rotate one of the samples so as to lift the spin orthogo¬ 
nality as shown in Fig. (3) right panel, which allows us to contrast the presence and absence of the orbital pseudo spin 
orthogonality by observing zero or finite current respectively in the junction. As shown in Fig. (3), the top and bottom 
surfaces with opposite orbital pseudo spin (r^,) = ±1 (denoted by the blue and red shades) have opposite spin textures 
even in the out of plane component. A rotation about x-axis by tt lifts the orbital pseudospin orthogonality but still 
the spin orthogonality persists. This is where the 6 -fold symmetry in the out-of-plane spin texture is exploited and 
it can be seen that a rotation of the second sample by 7 r /3 about the crystal growth axis, leads to a configuration 
where the states of the two surfaces with same momentum have identical spin components out-of-plane even though 
the in-plane components continue to be opposite. This is enough to create a finite overlap between the states in 
the spin sector and lift the spin orthogonality. In this situation, increasing the doping increases the warping effect 
leading to an increase in the overlap of the spin polarization of the two surfaces. Hence the tunneling current grows 
as a function of doping (see Fig. (4)). The fact that the finite tunneling current comes only from the non-zero overlap 
of the out-of-plane spin components can be verified by noting that the overlap of the spin components at momenta 
k* = n 7 r /3 is exactly proportional to the total current as shown in Fig. 4. For obtaining a geometric interpretation of 
the effect, maximally overlapping spins of the two surfaces state at momenta corresponding to states with maximal 
out of plane polarization (which contribute maximally to the tunneling current) is schematically shown on the Bloch 
sphere as a function of doping in inset of Fig. (4). 

So, we have explicitly identified that for 6 > = 0 , tt surfaces the orbital pseudospin (r) degree of freedom constitutes 
a half-metallic degree of freedom owing to complete polarization which is independent of momentum. We further 
showed that how the effect of warping can be exploited to demonstrate this half-metallic behavior via spin-valve-like 
effect, but with the orbital pseudospin. For arbitrary surfaces ^ 7 ^ 0, tt the half-metallic degree of freedom is degree 
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FIG. 4. The current, G;, (normalized by the current for \zi^\^ = 1 Vk) for Junction-B and \z\^* (^k* = 0) is plotted against the 

average chemical potential /i. The plot shows that they are linearly proportional to each other. The inset shows the behavior 
of the spins on the = 0 line for the two surfaces on the Bloch sphere for three different fi. 



FIG. 5. Left: The two spectrums, one conical and one ellitpical are shown on top of each other and the dotted line shows the 
locus of the points with equal energy and momentum. The subplots show the Fermi surfaces at four different energies, which 
show the presence of two or four, or the absence of Fermi surface intersections. Center: The normalized current (/) is plotted 
as a function of 0 for four different doping. To lift the no-match condition at higher negative dopings of the circular Fermi 
surface, one needs to deviate more from ^ = 0 as eccentricity of the elliptical Fermi surface increases with 0 till 7r/2 and then 
goes down. Right :The asymmetry in the current about ^ = 7r/2 is understood by looking at the contribution of the different 
sectors to the current. 


of freedom and similar calculations can be done to demonstrate it. Note that, since for the ^ = 0 , 7 r surfaces, the 
orbital pseudospin degree of freedom is completely frozen and has no dependence on the momentum (T^(k) = ±1), 
twisting the two samples in the junction about the crystal growth axis does not affect the orthogonality arising from 
the half-metallic behavior. On the other hand, for ^ 7 ^ 0 , 7 r surfaces, though the degree is frozen (T^(k) = ± 1 ) 
leading to the half-metallic behavior but in junctions of the form Junction-B the tunneling current could be affected 
by the anisotropy of the Fermi surface. Since the Fermi surfaces would be elliptical for 6 > 7 ^ 0, tt surfaces, twisting the 
two surfaces in the junction relative to each other can lead energy and momentum being simultaneously conserved 
only at discrete points leading to a different value of the tunneling current. For junctions of surfaces corresponding 
to different ^s the different shapes of Fermi surfaces (3) leads to further features in the tunneling current which we 
discuss in the next section. 


V. TUNNELING CURRENT BETWEEN ARBITRARY SURFACES 

Tunneling between the two opposite TI surfaces corresponding to the same 0 can always be suppressed due to 
the orthogonality in the half-metallic sector (T^) as was shown in Sec. IV. However, tunneling between two surfaces 
corresponding to different 0 can also be suppressed for an entirely different reason, namely the failure to conserve 
energy and momentum simultaneously, which we analyse in this section. 
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The Fermi surface of the surface states is actually an ellipse whose eccentricity depends on the angle 0 (3). Con¬ 
serving energy and momentum simultaneously would mean that the Fermi surfaces of the two TIs should intersect in 
the momentum space within the bias window. Failure to meet this condition leads to a zero current and we call it a 
no-match condition. The presence or absence of this condition is dictated by the angles the two surfaces make with 
their crystal growth axes and their respective dopings. For simplicity we take one of the surfaces of the junction to 
be the bottom surface perpendicular to the crystal growth axis with a circular Fermi surface doped by an amount D 
with repect to the second surface. If the second surface of the junction is allowed to be arbitrary (characterized by 0), 
then for a given energy the points of intersection of the Fermi surfaces denoted by k* = {k*,0k*) are parametrized 
by the doping fraction / {= D/E) and 0 as 


{ 7r/2,37r/2 if^ = -l 

0, TT if c = 1 (24) 

ei„ei, + Tr,0i^,ei^ + Tr if|ei<i, 

where ^ = 2 [v 3 /(l - /)]2 - [{vl + vl)/{vl - v|)]; 0%^ = cos-i(y(l + /)/2) and 01^ = cos-^{-^J(l + f)/2), with 
= {E — D)/hv\\. For |^| > 1, there exists no such intersections. 

This is illustrated in the left panel of Fig. 5 by overlaying the two spectrums and showing the presence or absence 
of intersections at different energies. The energy spectrum in Eq.(3) shows that for the elliptical Fermi surfaces, the 
semi-minor axes is equal to the radius of the circular Fermi surface which means they always touch at two points 

for D = 0, however when one of the samples is doped, the Fermi surfaces intersect only above a threshold energy. 
Using the locations of Fermi surface intersections from Eq.(24) and the overlap of the spinors at these intersections 
from Eq.(21), the current expectation is evaluated and normalized by a value that would have been present had 
the two spinors been parallel in the SU(4) space to define a dimensionless quantity (/) plotted in the center panel 
of Eig.5. The dimensionless current / is plotted as a function of 0 for different dopings of the first sample (with a 
circular Eermi surface) keeping the center of the tiny bias window fixed at a positive chemical potential well away 
from the neutrality point. The current is always zero for D > 0, due to the aforementioned no-match condition, as 
in this case, the radius of the circle of the Eermi surface for a given energy E shrinks by D/hv\\ and there are no 
intersections with the elliptical Eermi surface of the other surface. However for D < 0, the circular Eermi surface’s 
radius increases by \D\/hv\\, and as one deviates from ^ = 0 the eccentricity of the Eermi surface increases leading to 
intersections. The eccentricity of the Eermi surface being symmetric about 0 = tt/2 increases with 0 till 7r/2, so for 
higher negative dopings of the circular Eermi surface, one needs to go to higher deviations from 6> = 0 or tt to lift the 
no-match condition and get a finite current. This is observed the fact that range of 0 around 7r/2 over which there 
is finite current, shrinks for higher negative dopings. Also, note that for D < 0 the tunneling current is aymmetric 
about 0 = 7r/2. This can be understood from looking at the contributions of individual terms to the generalized 
tunnel magnetoresistance separately as shown in the right panel of Eig.5. 


VI. ROBUSTNESS OF THE HALF-METAL TO SURFACE POTENTIALS 

The half-metallic behavior of the surface state of 3D TI and its diagnostics via the orbital orthogonalities solely rely 
upon the fact that the surface states on the two opposite planar surfaces of a TI sample are orthogonal eigenstates of 
the Tq operator. However, a realistic scenario would naturally involve influence of surface potentials^^’^^ which need 
not keep this orthogonality intact. The surface potentials considered preserve time reversal symmetry(TRS), else the 
spectrum of the surface states itself can become gapped due to its influence and that will not be of interest for our 
study. These potentials not only change the nature of the surface states leading to a change in expectation values of 
the spin and the orbital pseudospin but they can also shift the position of the Dirac point in energy leading to doping 
effects^^. We systematically study the influence of these potentials on the pseudospin-valve effect for the 6> = 0,7r 
surfaces. 

The various classes of potentials allowed by TRS can be listed as I4 , (A^-T^) and {A.s'So)Ty which for the 0 = 0,7r 
surfaces look like I4 , (I2 •^) and (A^r • cr) where A^- is a vector that lives only in the x-z plane in the r-space 

and Act lives in the full 3D cr-space. It is important to note that the spin and orbital pseudospin degree of freedom are 
always decoupled for the 0 = 0,7r surfaces even in the presence of the surface potentials, so one can analyze their effect 
on the spin and orbital pseudospin sector separately. We first consider a surface potential of the form A S{z) I4 where 
the surface is taken to be at the z = 0 plane. This potential results in change of the orbital pseudospin texture of the 
surface states, but keeps the spin texture intact. As a result, the orbital orthogonality gets lifted. However we find 
the overlap in the orbital pseudospin sector ((1 + (t)i • (t) 2 )/ 2 ) is proportional to — Aq)‘^/{ m^ -|-Aq)^ where 

A = {Ao/m){2vz) and m is the bulk band gap of the TI. Note that the overlap tends to zero both in the m/Ao ^ 1 
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TABLE L The table shows the effect of the different surface potentials on the orbital (r) orthogonalities on both Junctions-A 
and Junction-B for the ^ = 0, tt surfaces. The column corresponding to 0 refers to case of no surface potentials. A /(X) denotes 
the presence (absence) of the orthogonality. 


and m/Ao ^ 1 limit and has significant contribution only in window around m/Ao ~ 1 / 2 . Hence for a junction 
with weak surface potential strength (A) will not effect the orbital orthogonalitys significantly. Next we consider the 
potential A S{z) l 2 0r^. Such a potential leaves both the spin and orbital textures of the surface states intact. Hence 
it does not effect the orbital orthogonality. Potentials of the form (g) I 2 change the orbital pseudospin texture of 
the surface states and hence it potentially can lift the orbital orthogonality. However, we find that the effect of this 
surface potential is the same on both the surfaces for Junction-A as far as the orbital pseudospin is concerned and 
hence for a symmetric junction (same strength of surface potential on both the surfaces) such a potential should not 
have any effect. Also, though this potential does not affect the spin texture of the surface states but it does lead to 
lifting of the spin orthogonality owing to shift in the Dirac spectrum in opposite direction on the two surfaces. Finally 
we consider potentials of the form A^r •cr( 8 )T^ which keep the orbital pseudospin texture intact, however they do affect 
the spin texture. Hence they preserve the orbital orthogonality expected to be observed in Junction-A, however in 
Junction-B since the orbital orthogonality is anyway absent, and these potentials destroy the spin orthogonality too 
and hence one can expect a finite conductance even in the absence of either an in-plane magnetic field or a doping to 
the warped regime of the spectrum. To summarize, we find that among all TRS preserving surface potentials, only 
(g) I 2 and I 4 affect the orbital pseudospin texture of the surface states. But even for these two potential we pointed 
out situations where the orbital orthogonality will not get destroyed. A summary of effect of these potentials on the 
orthogonalitys can be found in Table. I. 


VII. CONCLUSION 


We show that the surface states of 3D TI exhibit complete polarization in an emergent pseudospin (SU(2)) degree 
of freedom and hence behave like a half-metal. We find that this half-metallic behavior is robust against most of the 
surface potentials that preserve TRS. Since this half-metallic behavior leaves its signatures in transport properties 
between two surface states, transport probes to reveal this are also proposed. We showed that the half-metallic 
behavior can lead to a pseudospin-valve-like effect in the tunneling current between two TI surfaces, which can be 
used as a diagnostic for the half-metallic behavior. Our study indicates that weak disorder (which could get induced 
due to the roughness in surface potentials) can not only affect transport due to scattering of electron due to its 
coupling with the spin degree of freedom of the disorder potential (which is the one usually considered for study^^) 
but also due to its coupling with the frozen pseudospin degree of freedom of the disorder potential. Presence of such 
disorder may not be seen in the traditional methods used to probe surface states which couple only to the spin texture 
of the surface states such as the spin polarized ARPES and STM. Hence junctions of the type proposed in this letter 
can play a crucial role in understanding these pseudospin degree of freedom. It should be further noted that all the 
results obtained for the ^ = 0 , tt surfaces can be extended in a straightforward manner to arbitrary surfaces using 
Eq.(21) or Eq.(22) with a ^ So and r ^ Tq. For arbitary surfaces, the effect of the elliptic shape of the Fermi surface 
was also analyzed. 

To conclude, some remarks are in order about the applications of our results to topological insulator materials 
belonging to classes different from Bi 2 Se 3 . A notable example is the topological crystalline insulator in the SnTe 
material class^^’^^ which posses even number of Dirac cones on the surfaces. The surface states in such materials are 
not only spin-polarized but also have another degree of freedom namely the “mirror eigenvalue” whose topological 
robustness arises from the reflection symmetries of the crystal. We expect that the interplay of these different degrees 
of freedom can lead to further rich physics in the context of mesoscopic transport which is left here as a subject of 
future research. 
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